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Recent progress in the factorization of exclusive decays of heavy quarkonia into long-distance and short-distance 
contributions is presented. The role of higher Fock states and colour-octet contributions is outlined. 



1. INTRODUCTION 

Exclusive reactions at large momentum trans- 
fer Q can be calculated in perturbative QCD 
(pQCD) owing to a factorization theorem [JO that 
separates the short-distance physics of the par- 
tonic subreactions at the scale Q from the longer- 
distance physics associated with the binding of 
the partons inside the hadrons. The full ampli- 
tude is given as a sum of terms, where each term 
factors into two parts, a hard-scattering ampli- 
tude Th, calculable in pQCD, and wave functions 
tp(xi, kj_i) for each hadron H: 



(i.i) 



The separate dependence on the factorization 
scale fiF cancels, to the order considered, in the 
full amplitude A4. 



The importance of the various terms in (1.1) 
depends on their scaling with 1/Q. For the con- 
tribution that has the weakest fall-off with Q 
(leading- twist contribution), the amplitude Th 
describes the scattering of clusters of collinear 
partons from the hadron and is given by valence- 
parton scatterings only. Hence the only non- 
perturbative input required are the distribution 
amplitudes (DAs) <Ph {xi , Q) for finding valence 
quarks in the hadron, each carrying some frac- 
tion Xi of the hadron's momentum. The DAs 
represent wave functions integrated over trans- 
verse momentum k_L up to a factorization scale 
fiF, of order Q, and obey the BL evolution equa- 
tions 0. Indeed, most calculations so far have 
been performed in this standard hard-scattering 
approach (sHSA) summarized by: 
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• leading (i.e. valence) Fock states only; 

• collinear approximation; 

• (mostly) lowest-order in a s (Q); 

• tailored for S*-wave hadrons (orbital angular 
momentum L = 0). Up to now results of the 
sHSA remained at a qualitative level, lacking a 
rigorous discussion of (in particular the power- 
like oc 1/Q) corrections. Power corrections can be 
classified as follows: (i) Corrections arising from 
the overlap of the soft wave functions; (ii) correc- 
tions associated with the transverse momentum 
of the partons inside the hadrons; (iii) corrections 
from higher Fock states. 

Recently, exclusive processes at large scales 
have attracted new interest: 

• New, precise data on the pion-photon transi- 
tion form factor allow a quantitative extraction 
of the (non-perturbative) pion DA . 

• A modified hard-scattering approach (mHSA) 
has been developed, which allows for the incor- 
poration of transverse degrees of freedom and 
gluonic radiative corrections in the form of a 
Sudakov factor ||. This does not only al- 
low for a dynamical setting of the renormal- 
ization, a s (fim), and factorization, DAi(/Xfj), 
scales. Now also the perturbative contribu- 
tion to exclusive observables can be calculated 
self-consistently, i.e. without major contributions 
from phase-space regions, where the virtualities 
of the internal partons become soft and pertur- 
bation theory is not applicable. 

• A new factorization approach for inclusive 
quarkonium decays (and production) has been 
developed Q. Colour-octet contributions asso- 
ciated with higher Fock components are found to 





Figure 1. Feynman graphs for the colour-singlet 
decay \cj — ► tttt (J = 0,2). 



be crucial. 

• The importance of higher Fock states also for 
exclusive quarkonium decays has been realized 
and new power-counting rules derived . 



XcJ 



7T7T in the mHSA 



The decay XcJ into a pair of pseudoscalar 
mesons is given by the Feynman diagrams in 
Fig. [j], describing the annihilation of the domi- 
nant Fock state of \cj into a pair of gluons. For 
J PC = J++ mesons, the valence Fock state is a 
3 Pj colour-singlet cc pair (in the spectroscopic 
notation 2S+1 Lj of a spin angular-momentum 
state). Decays of heavy quarkonia can be cal- 
culated in the HSA since the heavy-quark mass 
m c provides the hard scale Q = m c . The sHSA 
result has been known for a long time || 
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where cj results from a convolution of the pion 
DAs with the hard kernel of Fig. [j], is the pion 
decay constant, and R' P (0) the derivative of the 
coordinate wave function at the origin. 

Recently, the calculations of the decay widths 
in the mHSA have been completed ||. The 
colour-singlet decay amplitude within the mHSA 
is calculated from the amplitude 



M (1) (Xc./->^) 



v/3 



J dxdy J d 2 bi d 2 b 2 exp[-5] 



H>*(yM)T H j(x,yMM)*AxM) ■ (2.3) 

This is the convolution w.r.t. the momentum frac- 
tions x,y and transverse separation scales bi,b2 



of the two pions; Thj is the Fourier transform 
of the hard-scattering amplitude Thj with the 
transverse (k±) dependence retained, see Fig. |]. 
The Sudakov factor exp[— S] takes into account 
those gluonic radiative corrections not accounted 
for in the QCD evolution of the wave function. 
Finally, the XcJ & re treated non-relativistically: 
the bracket in (2.3) contains the colour-singlet 
and spin projection operators, and ^j 1 ^ are the 
(reduced) non-relativistic wave functions, which 
yield Rp(0) after integration over the relative cc 
momentum k. The factorization and renormal- 
ization scales depend on the virtualities of the 
respective adjacent fermion propagators, e.g. 



f4u 



max <j 4(1 - x)(l - y) m 2 ,-^,-| 



(2.4) 



The results of the improvement of the mHSA 
calculation over the sHSA one can be summa- 
rized as follows: first, (soft) end-point contribu- 
tions are severely suppressed. Indeed, the results 
are self-consistent in the sense that only about 
2% of the decay rate originate from phase-space 
regions where a s (ii}ii)cts{^R2) > 1/2. Secondly, 
the uncertainty of the result against variation of 
the input parameters of the calculation (shape 
parameter B-2 of the pion-wave function parame- 
ter, charm quark mass, and QCD scale Aqcd) is 
reduced in the mHSA compared to the sHSA but 
still sizeable, see Table |l|. However, even with an 
optimistic parameter choice the theoretical pre- 
diction falls short by a factor of 4, calling for an 
additional contribution. 

3. Inclusive decays of heavy quarkonia 

The non-relativistic quark-potential model is 
very successful in describing the static proper- 
ties of heavy quarkonia. A quarkonium \H( J PC )) 
is considered to be pure, obviously colour- 
singlet, \QQ( 2S+1 Lj)} state where the QQ-pair is 
bound by an instantaneous potential. The non- 
relativistic nature of the bound system implies 
the existence of three important scales (besides 
the heavy-quark mass m): (the modulus of) the 
heavy-quark three momentum ~ mv, the quark's 
kinetic energy ~ mv 2 , and the QCD scale Aqcd- 

Inclusive quarkonium decays into light hadrons 
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Table 1 

Decay widths of Xcj ^ tt + it~ for various choices 
of the parameters compared with data. 



B 2 




Aqcd 




+7T-) [keV] 




[GeV] 


[GeV] 


J = 


J = 2 


Modified HSA 





1.5 


0.2 


8.22 


0.41 


0.1 


1.5 


0.2 


12.13 


0.53 


-0.1 


1.5 


0.2 


5.61 


0.33 





1.8 


0.2 


2.54 


0.12 





1.35 


0.2 


15.3 


0.78 





1.5 


0.15 


4.34 


0.21 





1.5 


0.25 


13.1 


0.68 


Experiment 




PDGjM 




105 ± 30 


3.8 ±2.0 




BES@] 




62.3 ± 17.3 


3.04 ±0.73 



are accessible to perturbative QCD: for suffi- 
ciently large m, the annihilation time ~ l/m 
is much smaller than the time scales relevant to 
QQ binding. Hence the short-distance annihila- 
tion of the QQ pair into gluons and light quarks 
can be separated from the long-distance effects of 
QQ binding. The basis of this separation is pro- 
vided by an effective field theory, NRQCD, sup- 
plemented by perturbative factorization 0. The 
inclusive decay rate is given as a double expansion 
in a s (m) and v 2 , to leading order in Aqcd/to- 

The factorization makes use of the fact that 
there arc, in any Feynman diagram, some prop- 
agators that are off shell by an amount ~ m 2 
much larger than the scales governing QQ bind- 
ing. Hence, in the effective low-energy theory 
these may be contracted to a point. The re- 
mainder of the diagram involves only distances 
much larger than l/m. For the case of inclusive 
decays, these are related to expectation values 
of universal long-distance matrix elements (MEs) 
O n - The O n are well-defined 4-fermion operators 
of NRQCD of the form O n = K n x)(^ K'J)) . 
The inclusive decay widths thus take on the form 



T = J2fn(H\On\H) , 



(3.5) 



The i mpor tance of the various long-distance 
MEs in (3J3) is given by their scaling with the 
heavy-quark velocity v. The velocity-scaling (po- 
wer-counting) rules were first established based 
on consistency requirements of the NRQCD field 
equations |J. More rigorously, these follow from 
a multipole expansion in the wavelength of the 
emitted (multiple) gluon field A = i?#/A g . In 
NRQCD, A ~ kg/(mv) ~ v. This formulation 
has the advantage that it can be generalized to 
the situation where Aqcd is not much smaller 
than mv 2 (as assumed in NRQCD) giving rise to 
different scaling rules |lC| ]. 

The Z/ip is dominantly an S-wave state. Its 
Fock-state expansion starts as 



|JMi— )) = o(i)|cci( a si)> 

+ 0(v)\cc s ( 3 Pj)g)+0(v 2 ) 



(3.6) 



where the subscripts 1 and 8 denote colour-singlet 
and colour-octet states, respectively. Annihila- 
tion into light hadrons via the valence Fock state 
involves the ME (J/^ |Ci ( 3 5'i ) | J/-0) , which we 
take to scale as v° = 1 (dividing out an over- 
all v 3 scaling). In the potential model this cor- 
responds to |i?s(0)| 2 . The first corrections for 
S-wave quarkonia arise from higher-dimensional 
operators (suppressed by v 2 ) still involving the 
leading Fock component. These lead, for exam- 
ple, to differences between the J/^ and rj c wave 
functions. The colour-octet component of ( pT^ ) 
enters first at relative order v A . Nonetheless, its 
contribution can be important, namely when the 
short-distance coefficients of the singlet states are 
suppressed by powers of a s (m) and/or kinemati- 
cal factors. An example is 3/tp production at high 
transverse momentum |fll"| . 

The situation is different for % c j, as becomes 
evident from their Fock-state expansion 



\Xcj(J- 



-+)> = 0(1) 
+ 0(v) |cc 8 



\<Zi(?Pj)) 

: 3 s 1 ) g ) + o(v 2 ) 



(3.7) 



where the short-distance coefficients /„ are calcu- 
lable as perturbation series in a s (m). 



The v 2 factor of the El multipole suppression of 
the ccg state is compensated by the fact that a P- 
wave operator scales as v 2 relative to an S'-wave 
operator. Hence, contrary to the case of S'-wave 
decays, two 4-fermion operators contribute to the 
decay rate of P-wave states into light hadrons at 
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Figure 2. Factorization of XcJ — ► tttt- 



leading order in v 

r[ XcJ ^LH] = ^{ XcJ \O^Pj)\ XcJ ) 

+ ^(Xcj|08( 3 ^i)|Xcj)+0( W 2 r) .(3.8) 

The term involving the colour-singlet matrix ele- 
ment is the one familiar from the quark-potential 
model. Indeed, up to corrections of order v 2 , one 
has (x,jWPj)\ Xc j) = 9|i?^(0)| 2 /(2^). 

The decays of P-wave (and higher orbital-angu- 
lar-momentum states) probe components of the 
quarkonium wave function that involve dynamical 
gluons. However, the contribution to the inclusive 
annihilation rate from the higher Fock state |ccg) 
is parametrized through a single number, namely 
the expectation value of the octet operator be- 
tween the XcJ state, i.e. the colour-octet matrix 
element in (|Q|). 

4. Factorization of exclusive quarkonium 
decays 

In contrast to inclusive decays, there are now 
two possibilities to absorb long-distance contri- 
butions, into either quarkonium-specific or light - 
hadron-specific objects. In order to set up a fac- 
torization scheme for exclusive quarkonium de- 
cays one has to consider the virtualities occurring 
in an arbitrary Feynman diagram. Three regions 
have to be distinguished. Propagators can be off 




Figure 3. Some colour-octet contributions to 

XcJ -> 7T7T- 



shell by order m 2 , (m,v) 2 , or Aq CD = (Am) 2 . Cor- 
respondingly, such contributions have to be asso- 
ciated with the hard amplitude, the quarkonium, 
or the light hadrons, respectively. This then leads 
to a triple expansion in a s (m), v, and A. 

Another difference with inclusive hadronic de- 
cays is the fact that the factorization into short- 
distance and long-distance contributions holds at 
the amplitude level rather than at the level of 
decay rates (i.e. squared amplitudes). In this 
sense factorization of exclusive decays resembles 
more the case of inclusive electromagnetic decays, 
where the final-state particles are also colour sin- 
glets. 

Higher Fock components in light hadrons cer- 
tainly exist. Consider the decay XcJ — * kit, 
Fig. ||. The (multipole-based) Fock-state expan- 
sion of the pion starts as follows 

|7T(0-+)) = O(l)|qq 1 ( 1 5 ))+O(A)|qq 8 ( 1 P 1 )g) 
+ 0(A 2 ){|qq 8 ( 3 S 1 )g) + ...}+0(A 3 ). (4.9) 
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Since both octet terms give corrections to the 
leading one of the same order A 2 , two new (non- 
perturbative) light-cone DAs had to be specified 
in a calculation using higher Fock components of 
the pion. Clearly, not much predictive power can 
be expected. In fact, higher Fock contributions 
to the pion0 are power-suppressed (oc Aq CD /m 2 ) 
and hence become negligible for sufficiently large 
quark masses. 

In contrast, higher Fock components of heavy 
quarkonia are suppressed by powers of v only, 
or even scale as the leading Fock contribution. 
Moreover, owing to the non-relativistic expan- 
sion one does not need new DAs but only new 
'decay constants' (expectation values of NRQCD 
operators). Of course, there remains the prob- 
lem of colour conservation. A perturbative fac- 
torization can be set up for the case a s (/j) /it <C 1 
for /i ~ mv since this allows one to treat the 
constituent gluon(s) as the constituent (valence) 
quarks. Colour conservation is then achieved by 
coupling the constituent gluon (in all possible 
ways) to the hard process, see Fig. ||. 

The calculation of XcJ ~ > i n the collinear 
approximation leads, however, to singular results. 
This is not surprising as the Feynman diagrams 
of Fig. H must contain contributions, which, to 
leading order in a s , constitute the two higher 
Fock states qqg of the pion. The regularization 
of the propagators associated with the 'would-be' 
constituent gluon of the pion depends on the hi- 
erarchy of A and v 2 : the singularities are pro- 
tected by either relativistic correction oc mv 2 to 
the heavy-quark momentum or the finite trans- 
verse momentum k± oc Aqcd = Am of the (light) 
valence quarks of the pion. 

In f§ v 2 < A (< v) was assumed, in which 
case it is the non-zero value of k± taken into ac- 
count in the mHSA that acts as a regulator and, 
in turn, determines the size of the colour-octet 
contribution to the decay width. The only new 
non-perturbative parameter is then a single num- 
ber, namely the octet decay constant /( 8 ) enter- 



2 The situation may be different for decays into higher-L 
mesons. 



Table 2 

Results for the XcJ decay widths (in keV) into 
pions (/^ = 1.46 x 1CT 3 GeV 2 ; B 2 = 0) in com- 



parison with experimental data. 









PDG § 


BES || 


XcO - 


-> 1T + 1T 


45.4 


105 ± 47 


64 ±21 


Xc2 ~ 


■* TT + TT~ 


3.64 


3.8 ±2.0 


3.04 ±0.73 


XcO - 


->7T° TT° 


23.5 


43 ± 18 




Xc2 " 


■+TT°TT° 


1.93 


2.2 ±0.6 





Table 3 



Results for the \cj decay widths into K's and 
77's in comparison with experimental data (/^ = 
1.46 x 10~ 3 GeV 2 , Bf = 0). 





J = 


J = 2 




r[xcj -> 


K+K"] [keV] 


Bf = -0.176 


22.4 


1.68 


Bf = -0.100 


38.6 


2.89 


PDG 1 


99 ±49 


3.0 ±2.2 


BES g 


52 ± 17 


1.04 ±0.43 




nxcj 


— * r]T]] [keV] 


B\ =-0.036 


24.0 


1.91 


Bl = 0.000 


32.7 


2.66 


PDG l! 


35 ± 20 


1.6 ± 1.0 



ing the |ccg) wave function 

\xf]) = f / (8) J d Zl dz 2 <£f(z l7 z 2 , z 3 )S^ (4.10) 

where t = A/2 is the Gell-Mann colour matrix 
and a the colour of the gluon. With this single pa- 
rameter XcJ decays into all pairs of pseudoscalar 
mesons can be predicted. The value obtained 
by fitting to the four decays into pions (Table ||) 
corresponds to a very reasonable probability of 
the octet |ccg) state, P C c g ~ 0.5. Predictions for 
XcJ decays into K's and 77's also agree nicely with 
data, see Table ||. Predictions for T decays have 
been obtained as well (Table |J) . 

It should be emphasized that the colour-octet 
contribution from the |ccg) Fock state is numeri- 
cally important, in accordance with the fact that 



G 



Table 4 

Decay widths in eV for P-wave bottomonia for 
two values of (m b = 4.5 GeV; R' P (0) = 
0.7 GeV 5 / 2 ). 



/W [GcV z ] 




— ► 7T + 7T 


Xb2 


— ► 7T + 7T 


1.46 x 10" 3 




20.6 




1.69 


5 x 10" 3 




152 




13.8 


/( 8 ) [GeV"] 


Xbo 




Xb2 


-> Tr°7T 


1.46 x 10" 3 




10.5 




0.88 


5 x 10" 3 




78.2 




7.27 



torization approach for exclusive decays of heavy 
quarkonia. More work is needed to generalize the 
factorization to higher orders in a s (mg), v, and 
Aqcd/"t-- In particular it has to be shown that 
the various contributions in an arbitrary Feyn- 
man diagram can be disentangled in a systematic 
way. 
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